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diameter  0.557  /«n.  Th*»  spectrum  of  the  scattered  light  is  a  very  narrow  Lorentzian 
whose  half  width  at  half  height  is  related  to  the  diffusion  constant  of  the  particles. 
The  scattered  light  is  allowed  to  fall  on  a  photomultiplier  tube,  and  the  times  be¬ 
tween  successive  photocounts  are  determined.  The  distribution  of  these  time  intervals 
is  not  random  for  times  on  the  order  of,  or  less  than,  the  coherence  time  of  the 
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could  be  expected  for  random  events.  The  "bunching”  of  the  photocounts  (photons)  is 
related  to  the  spectral  shape  of  the  scattered  light.  An  experiment  that  employs  the 
technique  of  pulse  interval  analysis  to  analyze  the  scattered  light  spectrum  is 
described  in  detail.  The  results  of  this  experiment  yielded  a  spectral  linevidth  that 
is  caused  not  only  by  the  diffusive  motion,  but  also  by  the  convective  motion  of  the 
particles.  This  determination  of  the  linewidth  was  found  to  be  in  agreement  with 
previous  experimental  results . 
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DIGEST 


The  purpose  of  this  work  is  to  study  the  light  scattered  from  a  simple  statistical  medium 
that  consisted  of  a  monodispeised  ensemble  of  spherical,  Brownian  particles.  The  light  from  a 
helium-neon  laser  operating  at  a  wavelength  of  632.8  nm  is  scattered  from  a  dilute  water  suspension 
of  the  particles,  which  are  polystyrene  spheres  of  diameter  0.557  nm.  The  spectrum  of  the  scattered 
light  is  a  very  narrow  Lorentzian  whose  half  width  at  half  height  is  related  to  the  diffusion  constant 
of  the  particles.  The  scattered  light  is  allowed  to  fall  on  a  photomultiplier  tube,  and  the  times 
between  successive  photocounts  are  determined.  The  distribution  of  these  time  intervals  is  not 
random  for  times  on  the  order  of.  or  less  than,  the  coherence  time  of  the  light.  The  probability  of 
counting  two  pulses  under  such  conditions  is  greater  than  could  be  expected  for  random  events. 
This  “bunching”  of  the  photocounts  (photons)  is  related  to  the  spectral  shape  of  the  scattered  light. 
An  experiment  that  employs  the  technique  of  pulse  interval  analysis  to  analyze  the  scattered  light 
spectrum  is  described  in  detail.  The  results  of  this  experiment  yielded  a  spectral  linewidth  that  is 
caused  not  only  by  the  diffusive  motion,  but  also  by  the  convective  motion  of  the  particles.  This 
determination  of  the  linewidth  was  found  to  be  in  agreement  with  previous  experimental  results. 
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I.  INTRODUCTION. 


%grrsg^rss 


A  LIGHT  SCATTERING  TECHNIQUE  TOR  STUDYING  BROWNIAN  MOTION 


Brownian  motion  is  a  topic  of  considerable  importance  in  the  study  of  aerosol 
mechanics.1  Many  diverse  aerosol  phenomena  involve  Brownian  motion.  Among  these  are:  (1)  dep¬ 
osition  of  aerosol  particles  from  a  stationary  gas;  (2)  dif.r-ion  of  aerosols  in  laminar  flow; 
(3)  filtration;  (4)  deposition  of  aerosols  in  the  respiratory  system;  (5)  absorption  of  aerosols  by 
bubbling;  and  (6)  coagulation  of  aerosols.  Also,  the  T  -w-aian  motion  of  aerosols  can  be  used  as  a 
particle  sizing  method.1 


There  is  a  50-year  tradition  of  extracting  physical  data  from  light  scattering  experiments, 
advent  of  lasers  with  their  very  narrow  linewidths  and  high  coherence  has  revolutionized  nearly 
all  aspects  of  light  scattering.  The  technique  of  studykj  Brownian  motion  by  means  >f  light 
scattering  is  one  area  that  has  been  opened  up  by  the  development  of  the  laser.  Although  this 
technique  has  so  far  been  applii  d  only  to  the  motion  of  particles  in  liquid  media,  the  application  of 
such  a  technique  to  aerosols  appears  very  promising. 


The  purpose  of  this  \  is  tc  *tudy  the  light  scattered  fro,  a  a  simple  statistical  medium 
that  consists  of  a  monodispett';<‘-  .nsemHle  of  ..physical,  Brownian  particles.  The  light  from  a 
helium-neon  laser  operating  at  a  v. .  ^engiJ  of  •''32  •<  >■■<:■.  is  scattered  from  a  dilute  water  suspension 
of  the  particles,  which  are  polystyrene  spheres  c<  d-.a?  f  ter  0.557  fim  (st:  figure  1).  The  spectrum  of 
the  scattered  light  is  a  very  narrow  Lorentzi.  ?  -  hose  half  width  at  half  height  is  related  to  the 

diffusion  constant  of  the  particles. 


HE-NE  LASER  LENS 


CELL 


Figure  1.  Experimental  Arrangement 


*  Fuchs,  N.  A.  The  Mechanics  of  Aerosols.  The  Macmillan  Company,  New  York,  New  York.  1964. 
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The  scattered  light  is  allowed  to  fall  on  a  photomultiplier  tube,  and  the  tinier,  between 
successive  photocounts  are  determined.  The  distribution  of  these  time  intervals  is  not  random  for 
times  on  the  order  of,  or  less  than,  the  coherence  time  of  the  light.  The  probability  of  counting  two 
pulses  under  such  conditions  is  greater  than  could  be  expected  for  random  events.  This  '‘bunching” 
of  the  photocounts  (photons)  is  related  to  the  spectral  shape  of  the  scattered  light. 


In  fact,  the  bunching  phenomenon  will  be  used  to  determine  the  spectral  linewidth  of  the 
scattered  light.  Such  a  pulse  interval  technique  is  a  very  attractive  method  of  studying  these  spectra 
because  as  the  linewidth  narrows,  there  is  a  corresponding  decrease  in  the  required  time  resolution. 

II.  BROWNIAN  MOTION. 


The  phenomenon  that  we  now  know  as  Brownian  motion  was  observed  for  the  first  time 
in  1827  by  the  botanist  Robert  Brown.  He  observed  that  material  particles  dispersed  in  a  liquid 
execute  an  incessant,  irregular,  and  randomly  directed  motion.  Subsequent  studies  by  several 
investigators  shewed  that  there  was  no  apparent  external  cause  for  those , notions.  At  the  beginning 
of  this  century,  Einstein  and  Smoluchowski  were  able  to  show  that  £rc*'man  motion  is  a  direct 
manifestationof  thermal  agitation;2*3  that  is,  the  motions  of  the  particles  are  produced  by  fluctuations 
in  the  collisions  with  the  molecules  of  the  surrounding  fluid.  In  their  approach  Brownian  motion 
was  treated  as  a  problem  in  random  walks. 

In  the  random  walk  approach,  a  Brownian  particle  undergoes  a  sequence  of  displace¬ 
ments,  s j ,  s2 ,  *  •  *,  tj,  because  of  collisions  with  the  fluid  molecules.  These  displacements  are 
Markovian  in  nature  because  the  magnitude  and  direction  of  each  displacement  are  independent^  of 
all  the  preceding  ones.  However,  the  probability  that  the  displacement,  sh  liesj>etween  %  and  s,+ 
dt(  is  governed  by  some  distribution  function,  Pify).  The  total  displacement,  S,  after  a  number  of 
such  displacements  is  given  by 


i 

The  solution  to  the  problemjs  given  by  determining  the  probability  /’(AS', A/)  that  the  particle 
acquires  a  net  displacement  AS  time  A/. 

If  At  denotes  an  interval  of  time  long  enough  for  a  particle  to  suffer  a  large  number  of 
displacements  but  still  short  enough  for  the  net  mean  square  increment  <IASP>  to  be  small,  then  it 
can  be  shown  that2*4 

rt.rto .  '0  >  =  l?" '° 1 1,4  ^  (2) 


where 


AS  =  r- 


A 


2Einstein,  A.  investigations  on  the  Theory  of  the  Brownian  Movement,  ed  R.  Furth.  Dover,  New  York,  New  York. 
1956. 

3Chandrasekhar,  S.  Stochastic  Problems  in  Physics  and  Astronomy.  Rev.  Mod.  Phys.  iJ,  1  (1943). 

4Wang,  M.  C.,  and  Uhlenbeck,  G.  E.  On  the  Theory  of  the  Brownian  Motion  II,  Rev.  Med.  Phys.  /7, 323  (1945). 


since 

r0  =  position  vector  of  the  particle  at  time  t  =  tQ 

f  =  position  vector  of  the  particle  at  time  t 
and 

A  t  =  t-  t0 

The  constant  D  is  known  as  the  diffusion  constant  and  is  defined  for  a  spherical  particle  that  obeys 
Stokes  law  as 


where 

k  =  the  Boltzman  constant 

T -  temperature 

77  =  viscosity  of  the  liquid 


(3) 


d  -■  particle  diameter 

Substituting  equation  2  into  the  definition  of  the  mean  square  displacement,  we  get  that 


<lr-r0P)  =  (47rDAt)"3/2  Pi 


Also,  by  definition 

<lr-?0 12>  =  <(x-x0)2)  +  ((y-y0)2)  +  <(z-z0 )2)  (5) 

And  because  the  liquid  medium  is  homogeneous  and  isotropic 

<(x-x0)2>=  l/3<lr-r0l2>=  2D  At  (6) 

The  equation  for  the  distribution  of  displacements  (equation  2)  and  the  equation  for  the  mean 
square  displacement  (equation  6)  have  both  been  verified  experimentally  by  Perrin.2'3 

III.  LIGHT  SCATTERING. 

Consider  the  cell  (see  figure  1)  to  be  filled  with  a  suspension  containing  M  identical  light 
scattering  particles.  The  incident  laser  light  is  polarized  perpendicular  to  the  scattering  plane,  and 
the  light  scattered  at  an  angle  9  is  observed  at  a  distance  R  from  the  illuminated  cell. 

Each  of  the  particles  is  subjected  to  an  incident  field,  £0e,<*o,',-woO>  where  k0  is  the 
wave  vector  of  the  incident  laser  beam,  and  Wg  is  the  corresponding  angular  frequency.  If  the 


Therefore,  the  phase  factor  is  just 


07=  lr7- r0  lcost?0  l4  I-  1^- r0lcost?,l*J=  (£0“£|>‘  ^"?o) 


For  convenience  let  us  define  the  scattering  vector, 


*  =  *o 


For  particles  of  velocity,  v  “<C,  we  have  that 


l*0  1*  lArJ 


and  from  figure  2 


where 


Iff!  =  2 l*0l  sin  0/2  =-^2  sin  0/2 


«0  =  refractive  index  of  the  medium 
X  =  wavelength  of  the  laser  radiation  in  vacuo 

-X 

In  terms  of  the  scattering  vector,  K,  the  phase  factor  is  given  by 

♦y  (/)«*•  (14) 

From  tjiis  relationship  we  see  that  (r)  changes  onl^if  the  particle  has  a  displacement  component 
along  K.  Any  displacement  of  the  particle  normal  to  K  leaves  fyO)  unaltered.  The  phase  factor  is,  of 
course,  a  random  variable  because  it  depends  upon  7^. 

The  quantity  - 1  Eq  I  depends,  in  general,  on  the  size  and  shape  of  the  scattering  particle. 
If  the  particles  are  spheres  that  are  small  enough  to  meet  the  Rayleigh  criterion,  cWX,  and  if  (m-  1 ) 
<  1 , then5 


^-j^j — -  •£<) ,sin  ♦  <w“  U 


where 


l£0'  \r  =  magnitude  of  the  electric  field  due  to  Rayleigh  scattering 


m  =  index  of  refraction  of  the  particle  relative  to  the  medium 

¥  =  the  angle  between  the  direction  of  E0  and  the  direction  of  propagation,  kt,  of 
scattered  wave 


5Clark,  N.  A.,  Lunacek,  J.  H.,  and  Benedek,  G.  B.  A  Study  of  Brownian  Motion  Using  Light  Scattering.  Amer.  J. 
Phys.  J0,  S7S  I V  .*n). 


If  the  light  scattering  particles  are  too  large  to  satisfy  the  Rayleigh  criterion,  but  they  do 
satisfy  the  requirement  that  (m  -  1)<1  and  the  Rayleigh-Gans  criterion,  kd/2(m  -  1 K1 » then  each 
volume  element  of  the  particle  gives  Rayleigh  (dipole)  scattering  that  is  independent  of  the  other 
volume  elements.  The  waves  scattered  in  a  given  direction  by  all  volume  elements  interfere  because 
of  their  spatial  distribution,  and  the  total  scattered  field  for  a  single  particle  is  found  by  integration 
over  the  whole  scatterer.  Under  these  conditions  it  can  be  shown  that* 


\Eo  \r.g  =  \}iV  f^  fdv] 


where  I  Eq  \r^  is  the  magnitude  of  the  electric  field  caused  by  Rayleigh-Gans  scattering  and  for 
spherical  particles 


'btdViw7)' 


Jm0K\d/2) 


t  I  i 

Equation  IS  gives  l£0 1^  for  the  case  of  spheres. 


The  index  of  refraction  for  polystyrene  spheres  relative  to  water  is  1.20  (sodium  D  line),6 7 
so  that  for  particles  of  0.SS7  jam  diameter: 

(m-  1 ) «  0.20  (18) 


kdl2(m-  0*0.55  (19) 

Therefore,  it  may  be  concluded  that  these  particles  are  scatterers  of  the  Rayleigh-Gans  type. 

Since  E( t )  and  K t)  =  l£(r)  I1  depend  upon  the  random  phases  of  the  waves  scattered  by 
each  particle,  it  is  impossible  to  predict  the  detailed  temporal  evolution  of  these  quantities. 
However,  such  a  random  process  can  be  characterized  by  quantities  such  as  Q( t)),  the  average 
intensity,  <£*(/)£(/  +  r)>,  the  correlation  function  for  the  fie]d,  and  <J( t)I(t  +  r)>,  the  correlation 
function  for  the  intensity.5*8  If  the  temporal  fluctuations  in  Eft)  and  Ift)  are  stationary,  then  the 
corresponding  correlation  functions  depend  only  upon  the  time  increment,  r.  The  spectrum,  Q( co), 
of  such  a  random  process  is  related  to  its  correlation  function  by  the  Wiener-Khintchine  theorem,4  -8 


Q(w)=  1/2 it  I  G(1)(r)<f/wTdr 


where  G(1)(r)  is  the  correlation  function  (G*(t)G(t  +  r)>. 


6Van  De  Hulst,  H.  C.  Light  Scattering  by  Small  Particles.  John  Wiley  &  Sons,  New  York,  New  York.  1957. 

7  Ashley ,  L.  E.,  and  Cobb,  C.  M.  Single  Particle  Scattering  Functions  for  Latex  Spheres  in  Water.  J.  Opt.  Soc.  Amer. 
4#,  261  (1958). 

8Kom,  G.  A.,  and  Korn,  T.  A.  Mathematical  Handbook  for  Scientists  and  Engineers.  McGraw  Hill,  New  York,  New 
York.  1961. 
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which  is  a  Lorentzian  centered  at  to  =  w0  and  with  a  half  width  at  half  maximum  of 

bio-D  \K  lJ  (29) 

With  the  aid  of  equation  3  and  the  fact  that  X  =*  632.8  nm,  the  linewidth  for  spherical  particles  in 
water  at  20°C  becomes9 


6w  =  98^0/2  (30) 

where  d  -  diameter  of  the  spherical  particles  in  pm.  From  the  correlation  function  (equation  26)  it 
can  be  seen  that  when  r  substantially  exceeds  the  time 

tc  =  1  /D\K\3  =  V8w  (31) 

the  scattered  field  approaches  zero.  This  time,  rc,  is  referred  to  as  the  coherence  time.  The 
expression  for  the  coherence  time  (equation  31 )  can  be  validated  formally  from  the  definition1 0 


I  C<1>(r)(*dr 


Thus,  an  important  result,  that  the  spectrum  of  the  light  scattered  from  a  dilute 
suspension  of  particles  undergoing  translational  diffusion  is  Lorentzian  and  contains  information 
about  the  motion  of  the  particles  including  the  diffusion  Constant,  has  been  reached.  For  the  case  of 
polystyrene  spheres  with  diameters  on  the  order  of  10* 1  pm,  an  extremely  narrow  linewidth,  bv  * 
102  Hz,  is  predicted  (equation  30).  Thus,  an  analysis  method  is  required  that  will  resolve  spectra 
that  are  too  narrow  for  the  usual  spectroscopic  techniques,  which  are,  at  best,  able  to  resolve  lines 
of  width,  bv  a*  107  Hz.  Such  methods  are  considered  in  the  next  section. 

IV.  SPECTRUM  ANALYSE. 

As  was  pointed  out  in  section  III,  the  spectrum  of  the  light  scattered  from  a  suspension  of 
dilute  polystyrene  spheres  necesitates  a  method  of  spectrum  analysis  that  is  more  sensitive  by  a 
factor  of  *  105  over  the  usual  dispersive  spectroscopic  methods.  The  interaction  of  light  with  a 
photoelectric  detector  provides  the  basis  for  two  related  techniques  of  successfully  measuring  such 
spectra:  (1 )  by  examining  the  photoelectron  current  spectrum  produced  by  light  beating;  and  (2)  by 
using  the  statistical  properties  of  individual  photoelectrons. 

A.  Light  Beating. 

Light  beating  may  be  described  as  the  mixing  at  a  photocathode  of  two  narrow  band 
optical  components  to  produce  an  a.c.  component  in  the  photocurrent  at  the  optical  difference 
frequency.  The  existence  of  this  phenomenon  was  first  demonstrated  experimentally  by  Forrester, 
Gudmundsen,  and  Johnson  in  1955. 11  Distributed  optical  spectra  also  can  be  resolved  by  this 


9Cununins,  H.  Z.,  Knable,  N.,  and  Yeh,  Y.  Observation  of  Diffusion  Broadening  of  Rayleigh  Scattered  Light.  Phys. 
Rev.  Lett.  12, 150(1964). 

10Mandel,  L.,  and  Wolf,  E.  Coherence  Properties  of  Optical  Fields.  Rev.  Mod.  Phys.  37, 231  (1965). 

1 1  Forrester,  A.  T.,  Gudmundsen,  R.  A.,  and  Johnson,  P.  0.  Photoelectric  Mixing  of  Incoherent  Light.  Phys.  Rev.  99, 
1961  (1955). 
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method.  Forrester1 2  considered  a  distributed  spectral  line  as  a  series  of  “slices”  and  found  by  a 
simple  extension  of  the  two-component  case  the  beating  between  all  possible  pairs  of  slices.  For  a 
Lorentzian  intensity  distribution  he  showed  that  the  power  spectral  density  of  the  photocurrent  is 
also  a  Lorentzian  centered  at  w  =  0  and  with  a  doubled  half  width.  Such  a  photocurrent  spectrum 
may  be  Fourier  analyzed  with  an  electronic  spectrum  analyzer  to  obtain  the  pertinent  spectral 
information. 


B.  Photoeiectron  Statistics. 


The  statistical  distribution  of  individual  photoelectrons  emitted  from  a  photocathode  is 
intimately  related  to  the  incident  optical  spectrum.  For  incident  light  with  a  perfectly  constant 
intensity,  the  probability,  pfn,7V,  of  counting  n  photoelectrons  in  time  T(T<  1/5 1>)  is 


pin,  T)  =  in)nln\e~<n>  (33) 

which  is  a  Poisson  distribution  with  variance 

<(An)J>  =  in)  (34) 

For  monochromatic  light  scattered  by  a  dilute  suspension  of  scatterers,  the  scattered  field,  ECr.t),  is 
a  Gaussian  random  process1 3  and  for  T<  1  /5 v 

pin,  T)  =  [( 1  +  <n»(  1  + 1  /<«»"]  " 1  (35) 

with  variance 


<(Aw)a>  =  <n>  +  <n>2 


(36) 


These  relationships  (equations  33,  34,  35,  and  36)  illustrate  the  general  situation  in  which  the 
photocount  distribution  will  depart  from  the  Poisson  form,  and  the  variance  of  the  number  of 
counts  is  in  excess  of  that  given  by  Poisson  statistics.1 0  The  nonrandom  arrival  of  the 
photoelectrons  is  a  manifestation  of  the  characteristic  bunching  properties  of  photons.  Several 
techniques  by  which  the  optical  spectrum  may  be  deduced  from  an  analysis  of  photoeiectron 
statistics  have  been  reported.14*18  The  use  of  delayed  photoeiectron  coincidence  measurements  to 
determine  the  linewidth1 9  *2 1  will  be  considered  in  detail  below. 

1  Forrester,  A.  T.  Photoelectric  Mixing  as  a  Spectroscopic  Tool.  J.  Opt.  Soc.  Amer.  51 , 253  (1961), 

1 3Cummins,  H.  Z.,  and  Swinney,  H.  L.  Progress  in  Optics.  Vol  VIII.  ed  E.  Wolf.  John  Wiley  &  Sons,  New  York,  New 
York.  1970. 

14Jakeman,  E.,  and  Pike,  E.  R.  The  Intensity-Fluctuation  Distribution  of  Gaussian  Light.  J.  Phys.  A.Proc.  Phys. 
Soc.  (London)  /,  128  (1968). 

15Mandel,  L.,  and  Wolf,  E.  Comments  on  a  Paper  by  E.  Jakeman  and  E.  R.  Pike,  ‘The  Intensity-Fluctuation 
Distribution  of  Gaussian  Light.’  J.  Phys.  A,  Proc.  Phys.  Soc.  (London)  1, 625  (1968). 

1 6Jakeman,  E.,  Oliver,  C.  J.,  and  Pike,  E.  R.  A  Measurement  of  Optical  Linewidth  by  Photon-Counting  Statistics.  J. 
Phys.  A,  Proc.  Phys.  Soc.  (London)  1 , 406  (1968). 

1 7  Arrecchi,  F.  T.,  Beme,  A.,  and  Sons,  A.  Measurement  of  the  Time  Evolution  of  a  Radiation  Field  by  Joint 
Photocount  Distributions.  Phys.  Rev.  Lett.  1 7, 260  (1966). 

18Be'dard,  G.  N-Fold  Joint  Photon  Counting  Distributions  Associated  with  the  Photoelection  of  Gaussian  Light. 


Phys.  Rev.  161, 1304(1967). 

1 9Morgan,  B.  L.,  and  Mandel,  L.  Measurement  of  Photon  Bunching  in  a  Thermal  Light  Beam.  Phys.  Rev.  Lett.  16, 
1012(1966). 

20Phillips,  D.  T.,  Kleiman,  H.,  and  Davis,  S.  P.  Intensity-Correlation  Linewidth  Measurement.  Phys.  Rev.  153, 1 13 
(1967). 

21Pusey,  P.  N.,  and  Goldburg,  W.  I.  Temperature  Dependence  cf  Concentration  Fluctuations  in  a  Binary  Liquid 
Mixture  Using  a  Photon  Correlation  Method.  Appl.  Phys.  Lett.  13, 321  (1968). 
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The  emission  of  an  electron  from  an  illuminated  photosensitive  surface  is  fundamentally  a 
quantum  mechanical  process  and  has  been  so  described  by  Glauber.2  2  For  the  optical  fields  under 
consideration  here,  however,  the  semiclassical  treatment  of  Mandel  and  Wolf10  is  entirely 
adequate.1 3  This  treatment  considers  the  optical  field  at  the  photocathode  to  be  specified  by  a  real 
vector  function 


m= r* 

J- OO 


oj)e~,utdai 


(37) 


Since  £<*>  is  real,  V(f,  w)  =  V*(r,  to),  and  E <r>  may  be  replaced  by  the  complex  analytic  signal2  3 


V(f. 


w)e~iutdo) 


(38) 


The  probability  of  photoelectron  emission  from  a  detector  that  is  illuminated  by  the  field  (equation 
38)  is  given  by 


p,  (t)St  =  o(E*(t)  ■  E(t)  St  (39) 

where  E*(t)  •  Eft )  is  the  instantaneous  intensity,  and  a  is  the  quantum  efficiency.  The  dependence 
of  px  (t)6t  on?  has  been  suppressed  on  the  assumption  of  spatial  coherence  at  the  detector. 

The  joint  probability,  p2(tx  ,t2)  St,  Sty,  of  registering  two  photocounts  at  times  t,  and 
r2  within  6^  and  6r2  is10-1  2 2 

p2(tx,t2)Stx  St  2  -  o2  <[£*(/,  )-E(tx)J  [E*(t2)‘E(t2)])St1  St  2  (40) 

For  a  stationary  field,  the  average  depends  only  on  the  time  interval  between  pulses,  t  -  t2  -  tx 
(figure  3).  Therefore, 

p2(/j  ,tx  +  t)  Stx  St  =  n2Ca)(r)St1  St  (41) 


l2 


Figure  3.  Pulse  Interval  Analysis  Scheme 


22Glauber,  R.  J.  The  Quantum  Theory  of  Optical  Coherence,  Phys.  Rev.  130, 2529  (1963). 
23Bom,  M.,  and  Wolf,  E.  Principles  of  Optics.  3d  Ed.  Pergamon  Press,  Oxford,  England.  1966. 
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where  the  normalized  second-order  correlation  coefficient  is  given  by 

(<2)(r,  _  <i£»(r  1 )  •  E(tj  )1  [g»(/|  +  r)  •  £(/,  +  r)] )  (42) 

1  <eu)£u))2 

and  h  -  mean  count  rate  of  the  illuminated  detector.  For  Gaussian  random  fields,  the  normalized 
second-order  correlation  function,  C<2>(r)  is  related  to  C(1)(t)  by5-13 

C<2>  (r)  =  [  1  +  lC<1>(r)l2l  (43) 

where  CG>  (r)  is  now  the  normalized  version  of  equation  26, 

C^1  Hr)  =  e~D  1  r  ■’  i  r '  c-/w,r  (44) 

Equation  41  assumes  that  all  counts  registered  are  a  result  of  the  absorption  of  photons  at 
the  detector  surface.  In  practice,  the  dark  current  is  responsible  for  a  mean  residual  count  rate,  m, 
which  is  unrelated  to  the  incident  radiation.  Equation  41  must,  therefore,  be  modified  to  read19 

p2(t,,f i  +r)5fjST  =  {Oi2  +m2)  +  n2 lC(1)(r)l2l  5tl8r  (45) 

The  number  of  pulse  pairs,  N,  that  will  be  counted  with  a  time  increment  r  within  the  range 
Tj  <t<t2  is  given  by 

fTj  • 

JV(rj,T2)  =  /s  (n2+M2)(r2~rl)+n2  I  |c<i>(r)l2  dr  (46) 

L  J*' 

where  ts  =  the  total  sampling  time. 

The  first  term  inside  the  brackets  of  equation  46  gives  the  random  rate  of  pulse  pair 
production,  which  is  constant  so  long  as  (r2  -  Tj )  is  unchanged;  the  second  term  represents  the 
“excess”  counting  rate  produced  by  the  bunching  phenomenon.  The  bunching  term  will  contribute 
significantly  only  when  t£  rc,  because  I  C^Ht)  I  (equation  44)  is  small  for  values  of  r  that  are 
appreciably  larger  than  rc.  For  this  reason,  the  upper  limit  r2  in  the  integral  (equation  45)  may  be 
approximated  by  r2  =  °°  so  that 


J  lC(i)(r)l2  dr- 


=  l/2Z?ltflV2Z>l*|,|T>1 


Applying  this  result  to  equation  46  gives 

Af(Ti )  =  /, l(«2  +  m2)(r2  -  T1)  +  n2/25we"2SwlT,J]  (48) 

According  to  equation  48,  a  determination  of  the  number  of  two  pulse  counts  as  a 
function  of  r2  will  provide  a  method  of  measuring  the  half  width  of  the  spectrum,  6w. 
Furthermore,  because  the  time  resolution  requires  decreases  as  5o>  becomes  smaller,  this  method  of 
spectrum  analysis  becomes  very  attractive  for  narrow  spectral  lines  such  as  the  scattered  light  from 
particles  undergoing  translational  diffusion. 


V.  APPARATUS  AND  PROCEDURES. 

A  block  diagram  of  the  laboratory  experiment  is  shown  in  figure  1 .  This  apparatus  allows 
us  to  measure  the  half  width  at  half  height,  6w,  of  the  scattered  light  by  determining  the  number  of 
two  pulse  counts  as  a  function  of  rt  (equation  48).  The  experimental  equipment  consisted  of  (see 
also  appendix  A)  the  following. 

A.  Laser. 

A  helium-neon  laser  with  ar.  output  of  10" 3  watts  at  632.8  nm.  The  laser  beam  is 
polarized  perpendicular  to  the  scattering  plane. 

8.  Lens. 

A  lens  of  7.5-cm  focal  length  was  used  to  focus  the  laser  beam  in  the  sample  cell.  It  is 
desirable  to  do  this  because  the  number  of  coherence  areas  (see  paragraph  D)  in  a  given  solid  angle 
decreases  as  the  scattering  volume  decreases. 

C.  Sample. 

A  spectrophotometer  cell  of  dimensions  1.0  by  1.0  by  4.7  cm  was  used  to  hold  the 
sample.  The  sample  was  prepared  by  diluting  the  polystyrene  emulsion  with  distilled  water  until  no 
halo  effect  from  secondary  scattering  could  be  observed.  Because  0.557-/xm  spheres  are  strong 
forward  scatterers,7  there  was  no  difficulty  in  eliminating  the  secondary  scattering  while 
maintaining  a  strong  scattered  beam  at  an  angle  of  30°. 


D.  Aperture  Stops  A  and  B. 

The  purpose  of  the  aperture  stop  A  was  twofold:  (1)  to  prevent  any  stray  light  scattered 
from  the  edges  of  the  sample  cell  from  reaching  the  photocathode;  and  (2)  to  provide  a  narrow 
beam  of  the  scattered  light  with  a  well-defined  scattering  angle.  The  aperture  A  was  circular  with  a 
diameter  of  0.071 1  cm. 


The  apert  B  was  located  at  a  distance  of  45  cm  from  the  sample  cell.  The  size  of 

this  aperture  is  of  .iportance  because  the  detector  area  must  be  limited  to  that  area  over 

which  the  optical  field  is  spat' ally  coherent  (see  page  12)i.e.,  the  area  of  coherence.  Detector  areas 
larger  than  the  area  of  coherence  will  produce  pulses  that  will  be  correlated  with  each  other  in 
random  fashion.  Cummins  and  Swinney1 3  give  a  detailed  calculation  of  the  area  of  coherence*  for  a 
quasi-monocnromatic  source  that  shows  that  for  a  single  coherence  area,  the  detector  area  should  be 


A<\2in 


(49) 


where  ft  =  the  solid  angle  that  the  source  subtends  at  the  detector.  For  this  experiment  aperture  B 
was  circular  with  a  diameter  of  0.071 1  cm,  which  gives 


A(K2ISlTl  *0.53 


(50) 


•See  also  Mandel  and  Wolfe;10  Forrester,  Gudmundsen,  and  Johnson;1 1  and  Forrester.24 
24  Forrester,  A.  T.  On  Coherence  Properties  of  Light  Waves.  Amer.  J.  Phys.  24, 192  (1956). 
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r.f.  =  10 -i!ClC1RlRL  *  0.015  (51) 

where  RL  -  total  resistive  load  following  the  filter  section.  The  resistor,  /?2,  serves  to  isolate  the 
pulses  appearing  at  the  anode  from  the  filter  section,  whereas  the  capacitor,  C*,  isolates  the  high 
voltage  d.c.  from  the  output. 

The  output  pulses  from  the  PMT  that  appear  at  the  anode  are  integrated  on  the  stray 
capacitance,  C,  which  appears  at  the  output  (figure  5a).  Let  us  consider  the  interaction  of  an 
idealized  pulse  (figure  5b): 

/0  =  0,  0  >  t  >  T 

i0  =  i,  +/2,0  <t<T  (52) 

with  this  circuit.  The  resulting  voltage  output  is  given  by  (see  figure  5c), 


Vs(t)=l/Cl  i1dt  =  i2R=i0R(l-e-,/RC),0<t<T  (53) 

Jo 

and 

Vs(t)  =  v0e~,/RC,t>T  (54) 

The  magnitude  of  the  stray  capacitance  is  on  the  order  of  100  pf  so  that  RC  «  10” 6  second. 
Because  T/RC  »  10"6,  a  first-order  expansion  of  equation  52  may  be  used  to  find  the  maximum 
voltage  developed  by  the  pulse 

V0  =  Ge/C  »  5  X  JO"4  volts  (55) 

where  e  =  electron  charge.  Therefore,  the  output  pulses  of  the  PMT  will  have  a  duration  of 
approximately  1  Msec  and  a  maximum  voltage  of  about  0.5  (mV). 

F.  Amplifier. 

The  multichannel  analyzer  accepts  pulses  of  0  to  2  volts  with  a  duration  of  approximately 
1  Msec  at  the  direct  access  connector.  The  selection  of  R  =  10KJ2  stretches  the  pulses  to  the 
required  time  duration,  but  amplification  on  the  order  of  103  to  104  is  required  to  raise  the  voltage 
to  the  desired  level. 

A  pulse  amplifier  circuit  (figure  6)  was  constructed  on  a  model  RP  operational 
manifold.28  This  manifold  is  provided  with  plug-in  receptacles  that  accept  up  to  eight  operational 
amplifiers.  A  color-coded  jack  field,  which  permits  the  installation  and  interconnection  of  circuit 
elements,  and  a  regulated  d.c.  power  supply  are  an  integral  part  of  the  manifold.  Four  amplifier 
stages,  each  employing  a  model  EP85AU  amplifier,27  were  constructed  in  a  simple  inverting 
configuration.  The  unity  gain  crossover  point  for  these  amplifiers  is  about  2  MHz;  because  no 
provision  is  made  for  changing  the  internal  compensation  of  these  amplifiers,  amplification  by 
successive  stages  was  used  to  preserve  as  much  bandwidth  as  possible. 


28Philbrick  Model  RP  Operational  Manifold.  Philbrick  Researches,  Dedham,  Massachusetts.  1967. 
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Figure  5.  Pulse  Output  Circuitry 


ftb  =  lOOKO  G  =  (f?b/Ra>4~:i04 


Figure  6.  Pulse  Amplifier  Circuit 


The  performance  of  the  amplifier  circuit  was  checked  by  inputing  the  calibrated  square 
wave  from  an  oscilloscope.  The  square  wave  inputs  were  0.5  mV  and  0.2  mV  PTP  and  had  a 
frequency  of  !03  IIz.  The  amplified  square  wave  pulses  were  observed  on  the  oscilloscope,  and  the 
predicted  gain  of  104  was  verified.  Some  distoiiion  in  the  form  of  narrowing  and  rounding  of  the 
square  wave  was  noted.  Amplified  dark  current  pulses  from  the  PMT  were  also  observed  on  the 
oscilloscope  and  found  to  have  peak  voltages  of  just  under  2  volts. 

G.  Multichannel  Analyzer. 

The  multichannel  analyzer  (MCA)  has  127  channels  that  are  available  to  count  voltage 
pulses,  and  it  has  the  ability  to  store  the  counts  in  memory  for  later  retrieval.  In  this  experiment  the 
MCA  was  used  in  the  multiscaling  mode  in  which  each  channel  will  accept  counts  for  a  certain 
preselected  dwell  time  and  then  transfer  to  the  next  channel  until  all  127  channels  have  been 
sampled.  Figure  3  illustrates  the  way  in  which  the  time  interval  fj  associated  with  a  2-pulse  count  is 
determined.  It  is  easily  seen  that 


tx  =/r0  (56) 

where  rc  =  the  dwell  time.  For  the  0.557-Mtn  polystyrene  spheres,  the.  coherence  time  for  a 
scattering  angle  cf  30'  is  tc  27  ms.  Therefore,  it  was  necessary  to  employ  the  shortest  dwell  time 
available  on  the  MCA,  1  ms,  to  execute  the  experiment. 


For  the  observation  of  the  spectrum  of  the  scattered  light,  the  apparatus  was  assembled  as 
shown  in  figure  1.  To  facilitate  optical  alinement  of  the  apparatus,  the  lens,  sample  cell,  and 
apertures  A  and  B  were  mou  .ted  on  small  laboratory  jacks.  The  operational  amplifiers  were 
balanced  in  accordance  with  the  manufacturer's  instructions,28  and  the  balancing  procedure  was 
repeated  at  half-hour  intervals  during  the  course  of  file  experiment.  The  scattering  angle  was  set  by 
making  the  distance  from  the  center  of  sample  cell  to  the  PMT  equal  to  45  cm  (±1  cm),  and  the 
distance  from  the  PMT  to  the  center  of  the  direct  laser  beam  was  made  equal  to  27.5  cm  (±1  cm). 
The  concentration  of  the  polystyrene  spheres  was  then  adjusted  so  that  single  pulses  were  being 
counted  by  the  MCA.  The  output  of  the  PMT  was  scanned  with  the  MCA  a  total  of  35  times.  At  the 
end  of  each  scan,  which  required  approximately  1 27  ms,  a  record  of  the  data  was  created  by  reading 
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the  memory  of  the  MCA  out  through  a  teletypewriter.  A  similar  set  of  data  was  produced  for  dark 
current  pulses  by  blocking  aperture  B  and  recording  the  data  as  described  earlier. 

VI.  RESULTS  AND  CONCLUSIONS. 

The  data  for  the  scattered  light  and  the  dark  current  were  analyzed  according  to  the 
scheme  presented  in  Section  V.  The  results  of  the  analysis  are  presented  in  figures  7  and  8. 

For  tiie  scattered  tight  an  estimate  for  the  optical  linewidth,  5a>.  was  obtained  by  a 
least-squares  fit  of  the  experimental  data  to  the  function  (see  appendix  B). 

N  =  A+BettTi  (57) 

where 

a  =  -26w 

The  parameters  A  and  B  depend  on  the  sample  time,  linewidth,  detector  area,  and  dead  time  effects. 
The  usual  procedure  of  treating  A  and  B  as  independent  parameters  has  been  followed  in  obtaining 
this  fit.29  The  least-square  fit  gives 

6o>  =  102  Hz  (58) 


k  2  3  4  5  6  7  8  9  10  X11213  14  1516  17  18  19  202122232425 

Tj  (ms). 


Figure  7.  Spectral  Data 

29Jakemin,  E.s  Pike,  E.  R,,  and  Swain,  S.  Statistic^  accuracy  in  the  Digital  Autocorrelation  of  Photon  Counting 
Fluctuations.  J.  Phys.  A:  Gen.  Phys.  3, 155  (1970 
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Figure  8.  Dark  Current  Data 


which  yields  a  value  for  the  diffusion  constant  of 

D  »  3.23  X  10“  7 cm2 /sec  (59) 

The  linewidth  as  determined  experimentally  is  considerably  larger  than  predicted  by 
equation  30,  which  gives,  5u  «  37  Hz.  This  excess  broadening  is  attributed  to  the  presence  of 
nonuniform  convection  currents  by  Cummins,  Knabler,  and  Yeh,9  whose  data  for  0.557-/um  spheres 
are  in  substantial  agreement  with  the  results  found  herein.  Because  no  attempt  was  made  to  isolate 
the  sample  thermally  and  acoustically,  some  improvement  in  the  results  may  be  possible.  The  use  of 
particles  of  smaller  size  would  render  the  effect  of  convection  broadening  relatively  less  important, 
but  the  time  resolution  required  to  perform  the  experiment  would  be  increased. 

The  results  obtained  from  the  dark  current  pulses  (Figure  8)  are  in  agreement  with  the 
prediction  of  equation  48,  that  the  random  rate  of  pulse  pair  production  is  constant  so  long  as  (r2  - 
Tj )  is  unchanged.  Therefore,  the  excess  coincidence  found  for  the  scattered  light  may  be  attributed 
to  the  motion  of  the  scatterers,  not  to  any  unusual  characteristics  of  the  dark  current  pulses. 

In  conclusion,  then,  we  have  seen  how  the  diffusive  motion  of  scatters  in  dilute 
suspension  produces  a  well-defined  spectrum  in  the  scattered  light.  It  was  shown  that  the  technique 
of  pulse  interval  analysis  provides  a  method  of  recovering  information  about  the  motion  of  the 
scatterers  from  such  a  spectrum.  Finally,  an  experiment  to  measure  the  spectrum  of  the  scattered 
light  was  designed  and  executed.  The  results  of  this  experiment  yielded  a  spectral  linewidth  that 
was  caused  not  only  by  diffusive  motion,  but  also  by  the  convective  motion  of  the  scatterers.  This 
determination  of  the  linewidth  was  found  to  be  in  agreement  with  previous  experimental  results. 
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For  our  purposes,  the  main  interest  in  such  light-scattering  studies  stems  from  their 
potential  application  to  particle  sizing  problems.  From  equation  30  it  can  be  seen  that  for  a  fixed 
scattering  angle,  a  determination  of  the  optical  linewidth  enables  one  to  find  the  Stokes  diameter  of 
the  particles.  However,  many  problems-such  as  the  effects  of  polydispersity  on  the  spectral  width, 
and  the  convective  motion  of  the  particles-must  be  solved  before  this  technique  becomes  practical. 
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APPENDIX  A 

EQUIPMENT  AND  MATERIALS  LIST 

1 .  Laser.  Spectra  Physics  Model  132-01 .  Spectra-Physics  Inc.,  Mountain  View,  California. 

2.  Sample.  Polystyrene  Latex:  0.557  nm.  Bioproducts  Department,  The  Dow  Chemical 
Company,  Midland,  Michigan. 

3.  Sample  Cell.  Opticell  Spectrophotometer  Cell.  Optical  Cell  Company,  Inc.,  Brent¬ 
wood,  Maryland. 

4.  Photomultiplier.  RCA  6217.  Radio  Corporation  of  America,  Harrison,  New  Jersey. 

5.  Power  Supply.  Atomic  Power  Supply.  Model  316,  Positive. 

6.  Amplifier.  Philbrick  Model  RP  Operational  Manifold,  Operational  Amplifiers:  Model 
EP85AU.  Philbrick  Researches,  Dedham,  Massachusetts. 

7.  Multichannel  Analyzer.  Nuclear  Data  Model  ND-1 10.  Nuclear  Data,  Inc.,  Schaumburg, 

Illinois. 


8.  Oscilloscope.  Tektronix  Type  545.  Tektronix,  Inc.,  Beaverton,  Oregon. 

9.  Laboratory  Jack.  Little  Jack.  Precision  Scientific  Co.,  Chicago,  Illinois. 
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APPENDIX  B 


LEAST-SQUARES  FIT  OF  EXPERIMENTAL  DATA 
In  order  to  fit  a  least-squares  curve  to  the  experimental  data,  the  equation  (equation  48  in 

text) 

N(fi  )  =  /,[(« 2  +  m2Xr 2  ~  tj ) +  /i2/25we_2wlr»  'J 

is  rewritten  as 

(N -  A)  =  BeaTl  (IB) 

where  a  =  -  26a) ,  and  A  and  B  are  regarded  as  independent  parameters.  Equation  1 B  is  then 
transformed  to  the  linear  form 

ln(N-  A)  =  In  B  +  a.rl  (2B) 

in  order  to  carry  out  the  least-squares  analysis. 

Because  equation  2B  involves  three  parameters  (A,  B.  ».■),  the  best  fit  is  obtained  by 
making  successive  estimates  of  A  and  calculating  B  and  a  in  the  usual  manner.  The  Gauss  criterion* 
is  employed  to  determine  the  set  of  parameters  that  give  the  best  fit  in  the  least-squares  sense.  This 
criterion  states  that  the  best  fit  is  obtained  when  the  quantity 

t—TFm-  <3B> 


where 

d(  =  deviation  of  the  measured  y  values  from  the  computed  values 


r  =  number  of  observed  data  pairs 


m  =  number  of  free  parameters 
is  minimized. 

Figure  B  shows  $  plotted  as  a  function  of  the  estimated  values  of  A.  It  is  readily  seen  by 
inspection  that  the  minimum  of  £  occurs  when  A  -  1 .00.  The  corresponding  least-squares  fit  is 

N=  1.00+  12Se_0,203rJ  (4B) 

where  Tj  is  in  units  of  milliseconds. 


•Worthing,  A.  G.,  and  Geffner,  J.  Treatment  of  Experimental  Data.  John  Wiley  &  Son*,  New  York,  New  York. 
1950. 
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Figure  B.  Minimization  of  the  Gauss  Criterion 
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